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Abstract. We report on the observation ol Bragg scattering ol an ultracold Fermi 
gas of ^Li atoms at a dynamic optical potential. The momentum states produced in 
this way oscillate in the trap for time scales on the order of seconds, nearly unperturbed 
by collisions, which are absent for ultracold fermions due to the Pauli principle. In 
contrast, interactions in a mixture with ^^Rb atoms lead to rapid damping. The 
coherence of these states is demonstrated by Ramsey-type matter wave interferometry. 
The signal is improved using an echo pulse sequence, allowing us to observe coherence 
times longer than 100 /is. Finally we use Bragg spectroscopy to measure the in-situ 
momentum distribution of the ®Li cloud. Signatures for the degeneracy of the Fermi 
gas can be observed directly from the momentum distribution of the atoms inside the 
trap. 

PACS numbers: 67.85.Pq, 37.25. 4-k, 03.75.Ss, 42.25.Fx, 37.10. Jk 
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1. Introduction 

Bragg diffraction of cold atoms at propagating standing light waves is used in two distinct 
ways, as Bragg spectroscopy to measure the momentum distribution of ultracold gases [1] 
and as beamsplitter for coherent atom optics [2]. Bragg spectroscopy has developed into 
a powerful tool for measuring the dispersion relation and the dynamic structure factor of 
cold gases. It was used to observe the momentum distribution of trapped Bose-Einstein 
condensates (BEC) [1, 3, 4, 5, 6], the structure factor of molecular condensates [7] and to 
study signatures of vortices [8, 9]. Furthermore, the technique is discussed for measuring 
the dynamic structure factor in the presence of long-range particle correlations, such as 
superfluid pairing in a Fermi gas or fermionic condensation [10, 11, 12]. 

Because Bragg diffraction coherently couples two momentum states, it is frequently 
used as a beamsplitter for matter waves in interferometric experiments [13]. A BEC 
seems the optimal candidate for such an interferometer, because the macroscopically 
populated wavefunction yields a high interferometric contrast. Consequently, a number 
of interferometric experiments have been done using Bragg diffraction of a BEC, e.g. 
[14, 15]. However, in the case of a BEC, the two momentum states resulting from the 
Bragg diffraction interact with each other via s-wave collisions [16]. This interaction 
is prominently observed in time-of-fiight (TOF) absorption images as a halo forming 
between the diffraction peaks. The typical lifetime of coherent superposition states due 
to interatomic interactions therefore is limited to only a few tens of microseconds [17]. 
Such experiments are thus performed during TOF or in very weak traps. 

Superposition states in a Fermi gas, however, can live for a very long time as has 
been observed for potassium atoms in an optical lattice [18]. Since fermions do not 
interfere, the matter wave contrast is necessarily limited to single-particle interference, 
but the long coherence lifetime nevertheless makes fermions excellent candidates for 
interferometric experiments. 

In this paper we report on Bragg diffraction of ultracold fermions from a light 
grating realized by two slightly detuned counterpropagating laser beams. The diffraction 
was analyzed by studying Rabi oscillations between the two coupled momentum states 
for low intensities of the Bragg lasers. Increasing the intensity of the Bragg light 
eventually leads to Kapitza-Dirac scattering. We also studied the transition between 
Bragg and Kapitza-Dirac scattering. Here, the system is no longer well described as 
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a two-level system and a multilevel ansatz has to be used. Experimental results of 
scattering in the Kapitza-Dirac regime and their theoretical description are presented. 

Furthermore, we used Bragg spectroscopy to map the momentum distribution of 
the ^Li cloud. As there is no interaction between identical fermions in the ultracold 
regime, the momentum distribution directly reflects the thermal distribution given by 
the Fermi-Dirac statistics. With a temperature of T/Tp = 0.6 the onset of degeneracy 
can be found by comparing Maxwell-Boltzmann and Fermi-Dirac fits to the measured 
momentum distribution [19]. 

The diffracted momentum states oscillate in the trap for several seconds nearly 
unperturbed, which is only possible because of suppression of s-wave collisions due to 
the Pauli principle. Higher partial waves are frozen out at the ultralow temperatures 
realized in the experiment. The oscillation is eventually destroyed by dephasing due 
to the anharmonicity of the trap. Furthermore, the lifetime of oscillating atoms is 
dramatically reduced in the presence of ^''Rb atoms that collide with the ^Li atoms. 

We also report on Ramsey interferometry with ultracold fermions, based on Bragg 
diffraction of trapped ^Li atoms. A complete Ramsey spectrum can be taken in a 
single experimental cycle and is recorded by absorption imaging after TOF. We observe 
Ramsey fringes for holding times of more than 30 /is. For longer times, the fringes 
cannot be resolved by the imaging system. 

The limitation can be overcome by echo techniques similar as in NMR [20]. This 
method has already been applied in atomic beam experiments [13] as well as in 
interferometric experiments with BECs [15]. In our experiment a diffraction echo is 
generated by a two-photon Bragg pulse located halfway in time between the Ramsey 
pulses. Hence, we realize a motional wavepacket echo by optical means in a similar way 
as it is done in a Ramsey-Borde interferometer [21]. Our paper extends previous work 
presented in [22]. 

2. Experimental procedure 

To achieve simultaneous quantum degeneracy of ^Li and ^^Rb, we use a procedure 
detailed in previous papers [22, 23]. ^Li atoms provided by a Zeeman slower and 
^^Rb atoms ejected from dispensers [24] are simultaneously collected by a magneto- 
optical trap. They are subsequently transferred via several intermediate magnetic traps 
into a loffe-Pritchard type trap, where they are stored in their respective hyperfine 
states \F,mF) = |3/2,3/2) and |2,2). This (compressed) trap is characterized by the 
secular frequencies {ujx,ujy,ujz)/2'7i = (762, 762, 190) Hz for ^Li and the magnetic field 
offset 3.5 G. For ^^Rb the trap frequencies are y^87/6 times lower. The ^'''Rb cloud is 
selectively cooled by microwave-induced forced evaporation and serves as a cooling agent 
for the ^Li cloud, which adjusts its temperature to the ^^Rb cloud through interspecies 
thermalization. However, due to the small interspecies scattering length, a = —20 as 
[23, 25], the thermalization is slow. Therefore, at low temperatures, when the size 
of the evaporated ^^Rb cloud becomes small or in shallow traps, where the different 



CONTENTS 



4 



gravitational sag separates the two species, the clouds thermally decouple. For the 
experiments described below, we typically reach temperatures of below 1 /xK with 2.5-10^ 
^^Rb atoms and 1.5 /iK with 2-10^ ^Li atoms. The respective densities are 2-10^^ cm"^ for 
®^Rb and 4-10^^ cm~'^ for ^Li, which corresponds to the critical temperature = 0.7 yuK 
and the Fermi-temperature Tp = 2.4 /iK. 




Figure 1. Schematics of a Bragg diffraction experiment, (a) Dispersion parabola 
of (quasi-) free particles. Shown is the initial Gaussian momentum distribution of 
the atoms and the Bragg transition between the initial state 1 and the final state 
2 with a momentum shift of 2hq. (b) Typical pulse sequence of a Bragg diffraction 
experiment. For single-pulse experiments only the first Bragg pulse is applied. For 
a Ramsey experiment two 7r/2 pulses are used. For a diffraction echo experiment, a 
third pulse (shaded area) is centered between the Ramsey pulses. 

For most experiments (unless stated otherwise) we use a decompressed magnetic 
trap with trap frequencies {ujx.ujy.ujz) /^ti = (236, 180, 141) Hz. Decompression lowers 
the temperatures of the clouds, but increases their gravitational sags. ^Li is much 
lighter than '^^Rb, so that the two clouds separate in space, and the ^Li cloud thermally 
decouples from the ^^Rb cloud. Since ^Li cannot thermalize by itself, the adiabatic 
cooling upon decompression becomes anisotropic. The temperature measured along the 
^-axis, along which the Bragg diffraction is performed, is then T^ ~ 0.9 /iK. 

2.1. Bragg diffraction 

The ^Li atoms are Bragg-diffracted by means of two counterpropagating laser beams 
aligned along the z-axis of the loffe-Pritchard trap. The laser beams are frequency- 
shifted by acousto-optic modulators (AOM). One AOM is driven by a stable quartz 
oscillator at 100 MHz. The other AOM is driven by a voltage-controlled oscillator, 
which is phase-locked to the quartz oscillator by means of an electronic feedback loop. 
The frequencies of the two laser beams are set to differ by an amount 6 = uj2 — uji = 
2hq^/m = 2ti ■ 295 kHz, with g = 27r/A and A is the resonant wavelength of the D2 line 
of ^Li. The Bragg beams have intensities of Ii = I2 = 13.7... 132mW/cm^. The laser 
frequencies are tuned A 1/2-11 ^ 1 GHz red to the D2 line, so that the two-photon Rabi 
frequency reaches values between Q-^ = Sttc^TI /Huj^Al ~ 27r ■ 47... 450kHz. The pulse 
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duration is chosen to generate specific pulse areas $ = nRr/27r. In most experiments, 
before applying the Bragg pulse, the ^^Rb atoms are removed from the trap by means of 
a resonant light pulse. After the Bragg pulse we wait for a time before switching off 
the trapping field and recording absorption images revealing the momentum distribution 
of the '^^Rb and ^Li clouds. In some (Ramsey) experiments, we apply a second Bragg 
pulse separated from the first one by a time interval At, or even a third one at half time 
between the two others (see figure 1). 




-10 12 

\ iq) 

Figure 2. TOF absorption images taken after a Bragg scattering pulse followed by a 
4 ms waiting time (half a trap period for rephasing) and a 2 ms ballistic expansion. 
From (a) to (d) the Bragg pulse duration is varied such that the pulse area is 
$ = ^2Rr/27r = 0, 0.4, 0.7, 1. 



2.2. Imaging 

Figure 2(a-d) show typical ^Li absorption images taken after application of a Bragg 
pulse. For low intensity of the Bragg beams, only a narrow slice is cut out of the 
fermionic momentum distribution. The position of the slice along the z-axis depends 
on the detuning of the Bragg lasers from the two-photon recoil shift, A = 2hq^ /m — 5, 
the width is due to power broadening by the Rabi frequency ^r. 

3. Theoretical model 

3. 1 . Free particles 

We describe our experimental observations with the following model [26]. We assume 
that during a Bragg pulse only two discrete momentum states j = 0, 1 of an atom are 
coupled and that the impact of the trapping potential can be neglected. This is justified 
for sequences much shorter than the trap oscillation period, i.e. r. At ^ 2tx/(jJz. The 
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probability amplitudes of the two states are denoted by . They correspond to atoms 
with initial momentum hkz that are coupled to states with momentum hk'^ = h{kz + 2q). 
The time evolution of the amplitudes under the action of the Bragg light is given by the 
solution of the Schrodinger equation 



with the Hamiltonian 



When the Bragg light is turned off, the Hamiltonian simplifies to 

By concatenating time evolutions described by e ^^^^/^ and e ^^^^^/^ the phase evolution 
of individual atoms in a superposition of momentum states can be calculated for 
arbitrary sequences of pulses separated by times of free evolution, e.g. Ramsey type 
sequences. 

Initially the atoms are inhomogeneously distributed in momentum space. This 
distribution is governed by Fermi-Dirac statistics. In the temperature range of our 
experiments {T/Tp > 0.6) however, the distribution does not deviate much from 
Maxwell-Boltzmann statistics, such that the momentum distribution can be assumed as 

To obtain the momentum distribution of the atoms after an applied pulse sequence, we 
calculate the evolution of the amplitudes aj^k^, ^01 a variety of initial momenta and weight 
the final populations of the momentum states with the distribution function (f){kz)- The 
number of atoms in the zeroth and first order Bragg-diffracted modes is then 

N,{t) = J (l)ikz)\a^,kMfdk,. (5) 

The procedure neglects atomic interactions, which certainly is a good assumption for 
an ultracold Fermi gas [27]. 



3. 2. Trapped particles 

If trapped atoms are considered, the problem arises that the momentum eigenstates are 
simultaneously coupled by two interactions, the moving optical lattice and the harmonic 
trap. However, the situation is simplified, if a separation of scales is possible. In general, 
the duration of a pulse is very short, r ^ 271 /uj^. In contrast, the duration of a waiting 
time period At can be such that it is no more negligible compared to a trap oscillation 
period. For these time intervals, the action of the trapping potential must explicitly be 
taken into account. 

Because the trap couples the atomic momenta and positions, the initial spatial 
distribution of the atoms must now be considered. For simplicity we describe it 
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as a thermal Gaussian function, similar as done in Equation (4) for the momentum 
distribution, 

Starting from initial positions z and momenta A;^, a diffraction pulse transfers the recoil 
2q to a part of the atoms. If the pulse length is short enough, the atomic position is not 
altered. The atoms then follow the classical trajectories 

kz it) = kz cos LUzt — z sin Uzt , (7) 

m 

k'Jt) = (kz + 2q) cos Uzt z sin Uzt, 

m 

where the first expression holds for undiffracted atoms and the second for diffracted 
atoms. These values are used instead of kz and k'^ respectively, in the Hamiltonians (2) 
and (3). As the Bragg pulses are short compared to the oscillation period, r ^ In jujz-, 
the effect of the trapping potential can be neglected for the description of the Bragg 
diffraction. This means that the Hamiltonian Hr depends on ^^(t), but can be treated 
as time-independent during the short time intervals r. In contrast, the free propagation 
Hamiltonian i^At gets time-dependent if At is long. In this case, the phase evolution 
of the atoms in both coupled momentum states can be described by writing the time- 
evolution operator as 

■exp(-z/o*rft^^K^)) 

^y:^[-^gdt[^k'^{t)~b 

Since the amplitudes aj,2,fc^ now also depend on the initial atomic positions, the 
final populations of the momentum states must additionally be weighted by the initial 
spatial distribution. Therefore, Equation (5) for the expectation values of the diffracted 
and non-diffracted atom numbers now reads 

^A^) = \\ <\>{kz)i){z) |a,,,,fc,(t)|'rfM^. (9) 



4. Single pulse experiments 

4.1. Bragg diffraction 

We first study diffraction with a single Bragg pulse. Having applied a pulse of a given 
length, we observe an axial modulation of the momentum distribution. This is due 
to the fact that different velocity classes of the atomic clouds have different Doppler 
shifts with respect to the Bragg lasers and therefore are subject to different Rabi 
flopping frequencies. The initial momentum distribution thus leads to an inhomogeneous 
population of the two momentum states after the pulse. 

For small pulse areas $ only momenta around hkz = are diffracted. With 
increasing $ nonresonant velocity classes are also coupled by the Bragg beams as can be 
seen in figure 3(a) and (c). The interpretation is supported by a theoretical simulation 
shown in figure 3(b) and (d), which has been done by computing the evolution of the 
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Figure 3. (a) Integration of the TOF absorption images shown in figure 2 along 
the radial direction (perpendicular to q). (b) Simulation of the axial momentum 
distribution after a Bragg pulse with same pulse area as in (a), c) False (color map of 
measured momentum distributions with pulse areas ranging between $ = f2Rr/27r = 
0...2.8. The radially integrated absorption images appear as rows, (d) Calculation of 
the false color map with the experimental parameters used for the measurement (c). 



momentum distribution according to Equations (2) and (4). Note, that the width of 
the initial momentum distribution exceeds the power and Fourier broadening in this 
experiment. 

The total amount of undiffracted and diffracted atoms is obtained by integration 
of the momentum distributions in the zeroth and first-order Bragg-diffracted clouds, 
respectively. The measured difference between the number of atoms in the first and 
zeroth order normalized to the total atom number is shown in figure 4 (crosses) as a 
function of pulse duration. The observed damped Rabi oscillations are reproduced by 
theory doing the same integration for the simulated momentum distribution according 
to Equation (5), as can be seen in figure 4 as a (red) dashed line. Note that the damping 
is solely due to the initial inhomogeneous momentum distribution, decoherence was not 
included in the simulations. 



4-2. Kapitza-Dirac scattering 
For larger Rabi frequencies. 



nK»^. (10) 

m 



with (Tfc^ = \Jmk-QT/h^ being the width of the momentum distribution, the Doppler 
broadening is dominated by power broadening, meaning that all atoms can be diffracted 
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Figure 4. DifFerence between the measured number of atoms in the zero and first- 
order Bragg-diffracted clouds for variable duration of the Bragg pulse. The (black) dots 
are the experimental data, the (red) dashed line is a theoretical simulation according 
to Equation (5). 



into the first order simultaneously. 

At some point, however, the Rabi frequency becomes comparable to the energy 
difference between adjacent momentum states, and Kapitza-Dirac scattering sets in. 
According to [26] the scattering will stay two-state like as long as the Rabi frequency 
fulfills the condition 

f^R < - (V - 2gafc J . (11) 

In the Kapitza-Dirac regime, the large energy uncertainty, connected with the fast 
coupling rate, allows several momentum states to be coupled simultaneously. The 
resulting dynamics can be described by a simple extension of the model detailed in 
Section 3. In particular the Hamiltonian (2) is replaced by 



a(fc^-2g)^ 



2m 
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(12) 



The crossover from the Bragg-diffraction regime to Kapitza-Dirac scattering is a smooth 
transition. As we will see below, for intermediate Rabi frequencies (here 1 MHz > 
nR/27r > 100 kHz) the neighboring diffraction states, corresponding to momentum 
shifts of 4:hq and —2hq (second and minus first order), are scarcely populated. For 
increasingly higher Rabi frequencies, Vt-^/2Tx > 1 MHz, the scattering populates more 
and more diffraction orders. 

We performed Bragg-diffraction at the transition between the Bragg and Kapitza- 
Dirac regime by increasing the Bragg laser intensity beyond the value used in 
the experiments discussed in the previous sections. Figure 5(a) shows the atomic 
distribution as a function of $ after 2 ms TOF for such an experiment. A simulation 
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using Equation (12) is shown in figure 5(b). This simulation reproduces the experimental 
results very well and reveals a Rabi frequencv of VLvt = 2tt ■ 420 kHz. 
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Figure 5. (a) False color map of measured momentum distributions similar to figure 3. 
The Rabi frequency was chosen such that the system is approaching the Kapitza-Dirac 
regime. While the main diffraction is between the momentum states of Ohq and 2hq 
(zcroth and first order diffraction states), diffraction into the —2fiq and Ahq states 
also occur. The pulse area was changed by applying different pulse lengths, (b) 
Theoretical simulation of the momentum distribution in (a) using Equation (12). The 
Rabi frequency chosen for the simulation was JIr ~ 2tt ■ 420 kHz. 

In some cases it is desirable to diffract all atoms simultaneously into the first order. 
This is only feasible by combining the Bragg regime (11), f^R <^ ^ (4g^ — 2qakJ with 
large power broadening (10), ^ 2hqakJ'm. The resulting condition, 



k^T < 



(13) 



is satisfied for temperatures T <^ 7 /iK. The temperature reached in this experiment are 
within this regime. 



4-3. Bragg spectroscopy 

For a given frequency detuning 5 the Bragg lasers address a specific momentum class 
hkz of the atomic cloud [1], 



5 



2hq^ ^ 2hqk, 



(14) 
m m 

This implies a linear dependence between detuning and addressed momentum class. 
It is thus possible to probe the momentum distribution by measuring the number of 
Bragg-diffracted atoms function of the detuning. 
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Figure 6. Momentum distribution of the trapped ^Li cloud measured by Bragg 
spectroscopy. The blue (solid) line is a Gaussian fit to the data, revealing a temperature 
of 1.5 /xK. The red (dashed) line is a simulated Fermi-Dirac momentum distribution 
for 2 • 10^ atoms at 1.4 /xK. These two curves are nearly indistinguishable in the plot. 

Figure 6 shows a spectrum resulting from such a series of measurements, performed 
in the compressed trap. In this experiment, the pulse duration was chosen such that 
the pulse area was $ ~ 0.5. The frequency detuning of the two Bragg lasers was 
changed between 40 kHz and 550 kHz in 15 kHz steps, which is small compared to the 
Rabi frequency ~ 27r ■ 58 kHz. In figure 6 the amount of scattered atoms is plotted 
versus the resonant atomic momentum derived from Equation (14). The data show a 
distinctive Gaussian shape, corresponding to the temperature T = 1.50 /iK ± 0.15 /iK 
(blue, solid curve), in good agreement with TOF measurements. 

At low temperatures T <^ Tp, the momentum distribution is expected to deviate 
from a Gaussian profile and to adopt the shape of an inverted parabola. However, in our 
experiments the temperatures were not far below the Fermi temperature of T-p ~ 2.4 /xK. 
In this regime, the distributions, integrated in the two dimensions perpendicular to the 
2;-axis, have quite similar shapes. A fit based on a Fermi-Dirac distribution for 1.8 ■ 10^ 
atoms is also shown in figure 6 as red (dashed) curve. This fit yields a temperature of 
T = 1.39 /iK± 0.18 /iK. 

The two curves, for the Gaussian and the Fermi-Dirac distribution are nearly 
congruent. Although the shape of the momentum distribution thus provides no clear 
signal for degeneracy in this regime, the result suggests that the temperatures due to 
Gaussian fits are overestimated by 5. .10%. 

Bragg spectroscopy yields the same information as standard TOF absorption 
imaging. However, this technique reveals the momentum distribution inside the trap. 
Adiabatic cooling and magnetic field distortions during trap switch-off do not affect the 
signal, which is particularly important for light and fast atoms such as lithium. 
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Trap dynamics 

By Bragg diffraction the atoms acquire additional momentum and begin to oscillate 
inside the trap. This oscillation has been analyzed in the following experiment. A 
fraction of the atoms, stored in the compressed trap, is Bragg-scattered by a single 
TT-pulse. The atoms are then held inside the trap for variable holding times before 
being imaged (see figure 7). 

Several observations are made: If no ^^Rb is in the trap, the ^Li atoms oscillate 
for times on the order of seconds; however, the distribution of the diffracted atoms 
smears out. Also, the oscillation period of the diffracted atoms differs from the period 
of undiffracted atoms. The momentum distribution of these undiffracted atoms is 
depleted around /c^ = 0, because atoms in low momentum states are Bragg-scattered 
preferentially. Atoms in the zeroth order oscillate with an average amplitude given by 
the width of the thermal distribution of the original cloud, as can be seen in figure 7. 




Figure 7. Momentum distribution after a 7r-pulse and a variable waiting time 
in the trap. Diffracted and undiffracted atoms oscillate in a potential given in good 
approximation by Equation (15). The equation of motion for this oscillation has been 
computed and is shown as (white and magenta) curves. The oscillation period for the 
two amplitudes differ, as can be seen for large waiting times 3> l-KjiOz- The (red) 
lines indicate the area in which the atom numbers were counted for analysis in figure 8. 



The broadening of the momentum distribution as well as the different oscillation 
periods for the two diffraction orders are due to anharmonicities in the magnetic trap. 
The magnetic potential of the trap is a combination of the quadrupole field generated 
by a pair of coils in anti-Helmholtz configuration and the field generated by four wires 
parallel to the rotational axes of the coils. The numerically calculated field geometry is 
in very good approximation given by [28] 

5^(^) = 5o-a(^ + a)+/3aV/^ (15) 

with the gradient a, the curvature /? and a = a/[3. With the resulting potential, the 
oscillation of the atoms in the cloud can be numerically simulated and compared to 
the anharmonic motion observed in the experiment. From this we derive values of 
a = 135 G/cm ± 25 G/cm for the gradient and /3 = 1.402 ■ lO^G/cm^ ± 4 G/cm^ for the 
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curvature, leading to the trap frequency ujz = ^JfJ^BPz/^ = 2vr ■ 181.6 Hz. The numerical 
calculation of the magnetic field configuration confirms the measured gradient [28]. 
Simple harmonic fits to the oscillation result in the same oscillation frequency with an 
error better than 10~^. This implies that the effect of the anharmonicity on the atomic 
motion is less than 1%. 
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Figure 8. Rephasing of the momentum distribution in the trap in the absence (a) and 
in the presence (b) of ^^Rb. Plotted are the integrated numbers of Bragg-diffracted 
atoms found at a given time after an initial Bragg 7r-pulse in a restricted area of the 
momentum space. The fit is a combination of an error function with linearly increasing 
Gaussian width and an exponential damping of the peak height. This damping is best 
observed in (c), where the integrated number of atoms in each peak is plotted versus 
time in a logarithmic plot. Only small damping is observed in the absence of ^^Rb 
(black, dash-dotted curve), while the presence of *^Rb induces a fast damping (red, 
dashed curve). 



We also analyzed the data by counting the number of atoms in a certain area 
of momentum space. This area was chosen as a box completely encompassing the 
diffracted atomic cloud at their highest momentum as indicated by the (red) lines in 
figure 7. The result of this analysis is shown in figure 8. Here, the fraction of atoms 
in the evaluation area relative to the total number of atoms is plotted as a function of 
time. The peaks are separated by a trap period and represent a revival of the original 
momentum distribution, directly after the Bragg pulse. For the fit in figure 8 (red 
dashed curve) we assume a Gaussian cloud shape oscillating in and out of the detection 
volume. The convolution of the Gaussian and the rectangular shape of the detection 
volume is described by a combination of two error functions with opposite slopes shifted 
with respect to each other by the time it takes for the cloud to oscillate in and out of 
the detection volume. The steepness of the slopes is given by the width of the cloud. 
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and this feature is repeated every trap period. 

As one can see from figure 8(a), the width a of tlie cloud's Gaussian momentum 
distribution broadens with time. If we assume a hnear broadening, the fit to the data 
yields 



where ctq is the initial width. This broadening is the main reason for the reduction in 
amplitude and is due to the anharmonicity of the trap. The integrated number of atoms 
in each peak, on the other hand, decreases much slower. This can be seen in figure 8(c), 
where the integrated number of atoms per peak is plotted against time. In the absence 
of ^^Rb we find a damping time of 3.2 s corresponding to more than 500 oscillations |. 

The exponential decrease of the atom number is due to p-wave collisions between 
diffracted and undiffracted atoms, since s-wave collisions are forbidden for a spin- 
polarized Fermi gas. The damping time therefore corresponds to a collision rate of 
7coii6p = 0.3 Hz. From this the p-wave scattering cross section can be calculated 
[23, 29]. Here we assume a Gaussian density distribution for the undiffracted atoms 
and a harmonic oscillation of the diffracted atoms. The relative scattering velocity 
= 2hq/m in this system corresponds to a scattering energy of 5 fiK ■ k-^. We derive 
a p-wave scattering cross section of cXp = 1.4 ■ 10"^^ cm^, which agrees with similar 
measurements performed on [27] and calculations made for ^Li [30]. 

In the presence of a cloud of ^^Rb, we find a dramatic reduction of the damping 
time, as can be seen in figure 8(b) and (c). The ^Li atoms may collide with '^^Rb atoms if 
these are not initially removed from the trap. The collision time constant calculated from 
the heteronuclear s-wave scattering length [23] is 'JcMgsy — 107ms. The exponential 
damping in the fit of figure 8(b) and (c) resulted in a damping time of 113 ms, while 
the broadening is of the same order as in the case without ^^Rb. There is excellent 
agreement between the calculated collision rate and the observed damping. 

5. Atomic coherences 
5.1. Ramsey inter ferometry 

The observed revivals of the density distributions after long times indicate that the ^Li 
atoms are free from any kind of perturbation, which in turn suggests long decoherence 
times for quantum mechanical superpositions. The existence and the stability of such 
states can be tested by interferometric experiments. 

To this end, we apply a Ramsey sequence of two 7r/2 pulses, like in an atomic Mach- 
Zehnder (or Ramsey-Borde) type interferometer. The beamsplitters are Bragg pulses 
with a duration equivalent to a pulse area of $ = 0.25, and the interferometric paths 

X In the decompressed trap, we found damping times without ^^Rb atoms exceeding 13 s, which is on 
the order of the magnetic trap hfetime. However, in those traps the presence of a ^^Rb cloud has no 
impact on the ®Li because of gravitational sag. 




(16) 
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are closed by the kinetic energy released during time-of- flight. We proceed as follows. 
At time t = in the decompressed trap, a first 7r/2 pulse of duration r is applied. Then 
we wait for a variable amount of time At after which a second 7r/2 pulse is applied. The 
momentum distribution is recorded after the delay of half a trap cycle and 2 ms ballistic 
expansion. Typical time-of-fiight absorption images are shown in figure 9 (a) and (b). 
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Figure 9. (a,b) Absorption images taken after a Ramsey pulse sequence and 2 ms 
ballistic expansion. (c,d) Integrations of the absorption images perpendicular to k^. 
The (black) solid line is the experimental data, the (red) dashed line is a theoretical 
simulation using Equations (2), (3) and (4). The evolution time between the two 
Ramsey pulses is At ~ 6 /is in (a,c) and A< = 12 /is in (b,d). The temperature used 
for the simulation is T = 1.2 /iK and the Rabi frequency is il^ = 27r ■ 58 kHz. 



The initial momentum distribution of the atoms gives rise to an inhomogeneous 
Doppler shift, which detunes the atoms from Bragg resonance. This allows to observe 
all Ramsey fringes in a single shot. We do observe Ramsey fringes for time separations 
of the Ramsey pulses on the order of up to 30 fis. For longer times, the fringe spacing 
falls below the resolution limit of the imaging system. 

Using Equation (3) the phase difference for the two momentum states of an atom 
moving at velocity hkz/m accumulated after a time At is given by 

^ 2h{q + h)q 



-At-S. 



m 



(17) 



Adjacent fringes are produced by atoms whose velocities differ by an amount hAkz/m 
such that A(j) = 27r, yielding the condition 
hAk, 



m 



vr 
'qAt 
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After ballistic expansion time during a time ttof the fringes are separated in space by an 
amount irttof/q^t, which must exceed the Ax = 14 /im resolution of our optical system. 
From this follows the condition for the maximum allowable Ramsey pulse separation, 

AWx = -^ = 48/is. (19) 
Ax ■ q 

The realistic value is additionally subject to noise. 

To a certain extent the contrast can be enhanced by longer TOFs. However this also 
dilutes the cloud and reduces the quality of the absorption image. Lower temperatures 
may allow for longer TOFs, but as the temperature of the cloud is below the Fermi 
temperature already, the width of the momentum distribution cannot be reduced by 
much due to the Pauli exclusion principle. 

The maximum detectable evolution time could be slightly improved by a better 
imaging system, however, the general problem remains. It can be circumvented by an 
echo technique described in the following section. 



5.2. Diffraction echo interferometry 

Interferometry of particle ensembles is subject to decoherence and dephasing. 
Decoherence destroys quantum mechanical superpositions due to coupling to a reservoir. 
In contrast, dephasing results from divergent phase space trajectories for different 
particles and may be reversed if the trajectories can be inverted [31]. This is possible 
with echo techniques well known in NMR spectroscopy [20] . They have been successfully 
applied in cold atom experiments [13] and BEG interferometry [15]. In an echo 
experiment the free evolution time between two Ramsey pulses is split in two intervals 
by an additional vr-pulse inverting the phase evolution of atomic states. Diffraction echo 
interferometry is robust against dephasing. 

In our case the signature for diffraction echo is the almost complete revival of the 
zeroth order momentum state. Such an experiment is shown in figure 10. It is similar 
to the Ramsey experiment described above, but the waiting time between the pulses is 
now divided into two waiting times ti (20 /is in figure 10(a)) and T2, separated by an 
additional vr-pulse. After completion of a pulse sequence, the atoms are held in the trap 
for half a trap period, then the trap is turned off and the atoms are imaged after 2 ms 
ballistic expansion time. 

The Rabi frequency in this experiment is Q-^ = 27r ■ 220 kHz. This ensures Bragg 
resonance for all atoms according to Equation (10). However, for this Rabi frequency 
Kapitza-Dirac scattering already occurs, as can be seen in figure 10. For times T2 ^ t\ 
the zero-momentum revival is incomplete, it appears as a central peak (dip) in the zeroth 
(first) order Bragg-diffracted portion. This peak (dip) broadens until for = ti the 
revival is most pronounced. This effect is also seen in figure 10(b), where the difference 
between atom numbers in the zeroth and first diffracted state is shown. 

The revival is seen for times up to rtot = ti + r2=100/is. In a harmonic trapping 
potential momentum states are no eigenstates, which leads to a mixing of position and 
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Figure 10. (a) Integrated density of absorption images taken after a spin echo pulse 
sequence and 2nis TOF, for a waiting time ti = 20 /is and a variable waiting time T2. 
(b) Deviation between the total number of atoms in the zeroth and first exited state. 
The theoretical curves in red (dashed) are simulations based on Equation (8), using a 
Rabi frequency of fi^ = 27r • 220 kHz. 



momentum coordinates for long evolution times. This process is responsible for the 
decay of the diffraction echo rather than decoherence. The coherent theoretical model 
based on Equation (8) fully describes the experimental observation, which leads us to 
claim Ttot as a lower bound for the coherence time. 

Setting the waiting times to a multiple of the oscillation period ti = T2 = 27m /uz 
can compensate the effect of the trapping potential as has been shown in [32]. In our 
experimental setup, revivals of the diffraction echo could not be observed. This might 
be due to the large momentum spread of the initial cloud. 

6. Conclusion 

In conclusion, we presented Bragg scattering of fermions from a moving optical lattice. 
We were able to show the coherence between the zeroth and first order diffracted 
momentum states. Our work presents the first study of interferometry with fermionic 
atoms using separate interferometric arms and a variable evolution time. The coherence 
is also observed with a diffraction echo experiment, where a lower bound for the 
coherence time of ttot > 100 /is was found. A longer observation was hampered by 
the fact that the studied momentum states are no eigenstates of the harmonic trapping 
potential. The diffracted atoms, however, oscillate in the trapping potential for very long 
times with a single atom damping time of about 3.2 s. This leads us to the assumption of 
a similarly long coherence time, which is consistent with the observation made by Roati 
et al. [18]. Thus fermionic ^Li is a promising candidate for interferometric experiments 
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with long coherence times using schemes similar to [15]. 

Moreover, we successfully applied Bragg spectroscopy to fermions, measuring the 
momentum distribution of an ultracold degenerate Fermi gas. Although not deep in the 
degenerate regime (T/Tp = 0.6), the momentum distribution revealed signs of quantum 
degeneracy. For even lower temperatures it is conceivable to observe degeneracy from 
the shape of the distribution directly. In addition, Bragg spectroscopy allows for in-situ 
probing of specific momentum classes of trapped atoms. 

Future work will be devoted to studying the ®^Rb+^Li mixture in a crossed 
beam dipole trap, using Bragg scattering in the vicinity of one of the recently found 
heteronuclear Feshbach resonances [33]. Signatures of the interaction energy of the 
ensemble should be visible in the profile of a Bragg spectrum, similar to the frequency 
shift observed in Bragg spectroscopy on BECs [1]. Furthermore, Bragg spectroscopy 
permits to measure the atomic dispersion relation, which is particularly interesting in 
case of many-body correlations. Bragg spectroscopy should also be applicable to take 
excitation spectra of fermionic spin mixtures in the BEC-BCS crossover regime which 
could be realized in our experiment as well. 
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